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Abstract 
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1 Introduction 

Fractional fields have often been used to model irregular phenomena. 
The simplest one is the fractional Brownian motion introduced in [5] 
then developed in [T3]. More recently, many fractional processes have 
been studied, usually obtained by a stochastic integration of a deter- 
ministic kernel against a random measure (cf. among others [H], [T], 
[7], [9] and [2]). Many different simulations methods were discussed in 
the literature, but shot noise series seem to fit perfectly in that kind 
of problem. Generalized shot noise series where introduced for simu- 
lation in [TB], further developments were done in [T7] and [TB] and a 
general framework was developed in [5]. Moreover, a computer study 
of the convergence rate of LePage series to a-stable random variables 
has been done in a particular case in |6]. 
A shot noise series can be seen as: 

oo 



where (Vn) are i.i.d. random variables and (r„) the arrival times of a 
Poisson process. Usually, there is no question about the simulation of 
Vn- In this paper, we will contribute to this theory by looking at the 
convergence rate when Vn can not be simulated but only approximated. 
Moreover, we will show a way to save computer time. Indeed, in a shot 
noise series representation of a a-stable process, the first terms are 
the bigger so you have to minimize the error made in approximating 
Vn for small n. For large n, Tn^^" is small that it is not as useful 
to approximate Vn them with so much care. We developed a way of 
measuring the convergence rate towards the limiting process depending 
on the error made while computing each term. 

Afterwards, we will consider an application of this result to study 
a particular process. In network traffic modeling, properties like self- 
similarity, heavy tails and long-range dependance are often needed, 
see for example [2]. Moreover, empirical studies like [5] showed the 
importance of self-similarity and long-range dependance in that area. 

In [4], the authors introduced fractional Brownian motion local 
time fractional stable motion as a stochastic integration of a non- 
deterministic kernel against a random measure, which will be our main 
interest in the second part. Here we will call it Local Time Fractional 
Stable Motion (LTFSM). This process was defined as: 

/ / l{x,t){Lu')M{duj',dx), fori>0. 
Jw Jr 

In this expression, I is the local time of a fractional Brownian motion 
of Hurst parameter H defined on {W, !F' , P'). M is a symmetric alpha 
stable random measure (see [12] for more details) with control measure 
P' X Leb (Leb being the Lebesgue measure on M) . LTFSM is a-stable 
but also self-similar and its increments are long-range dependant. 

The first step towards understanding LTFSM is naturally to have a 
look at its sample paths. Unfortunately, the above expression does not 
directly give a way to get the sample paths. In the case of Brownian 
motion local time, this process can be seen as the limit of a discrete 
random walk with random rewards model. It is not completely satis- 
fying for a few reasons: first, it only works for H — 1/2, then, there 
is no control of the convergence speed rate towards the limit. That is 
where the tool we developed in the first part comes to play. In this 
paper, we will study how we can simulate this process by using the 
expression given in equation (5.3) in [¥J, which can be seen as a shot 
noise expansion. 

The next section is devoted to the shot noise theory results. In the 
following part we will see how LTFSM fits in the general frame we just 
developed, and we will be able to get a convergence rate in the case 
of our process. The last section will be devoted to a quick study of 
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our simulations with a comparison to the random walk with random 
rewards model that was introduced in |4j in the case H — 1/2. 



2 Shot Noise series 

In this section, we will show some results on shot noise series, using 
theorem 2.4 in [17]. 

Assumption 2.1. Ek will he the space of continuous Junctions defined 
on a compact subset K <zM., equipped with the uniform norm denoted 
\\-\\k- ^6 will denote by IHj^p the LP{K) norm. Let us consider the 
case h :R X Ek — > Ek where 

h{r,v) ^r-^l'^v. (1) 

h is a Borel measurable map. Let (r„)„>i be the arrival times of a 
Poisson process of rate 1 in M^, (14i)n>i be a sequence of i.i.d. sym- 
metrical random variables of distribution A. Let us assume (r„)„gN, 
(l^n)neN oltc independent. 

We will suppose that the assumptions of corollary 2.5 in [T7] are 
satisfied, i.e. : 

Assumption 2.2. The measure defined for all A G Be (Borelians of 
Ek) by 



F{A)^ / / lA\{^}{h{u,v))\{dv)du (2) 
Jo JEk 

is a Levy measure. 

Now, we will get some inspiration from the proof of [lOJ to prove 
the following proposition, considering the L'^ distance between the sum 
of the series, and the truncated sum. 



Proposition 2.3. Under assumptions \2.1\ and \2.^ let us consider the 
shot noise series 



y = ^/i(r„,K). 



n=l 

Let us denote Yj^ — X^^iLi ^(-TjuKi)- Take q > 2, we have to suppose 
that there exists Mq such that for all n, ¥\\\Vn\W] — ^-^i Then, 
for N > q/a — 1, we have 

E[\Y{t)-Y^m<^^^^^, 

where Aq depends only on q through Mq and Bq, and 

Tin - ^) 

" r(n) 

with lim„^oo Hn^q = 1. 
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Proof. Denote for < P 



n=N+l 



According to theorem 2.4 in |17j . limp_>+oo Rn,p exists. 

As h{Tn,Vn){t) = Tn^^"Vn{t) and because Vn{t) is symmetric, 
h{Tn,Vn){t) is also symmetric. Therefore, we can apply proposition 
2.3 in [TT], and we get 



E 



max \Rn nit)\'' 



<2E[|i?jv,p(i)|' 



We now are able to apply Khintchine inequality Let e„ be a sequence 
of i.i.d. Rademacher random variables, independent from everything 
else. Thus e„/i(r„, V^)(t) has the same distribution as h{Tn,Vn){t) 
since it is symmetrical. Khintchine inequality claims 



E 



^ e„/i(r„,K)(i) 



n=N+l 



(Mr„,y„)(t))„eN 



1/9 



< 



/ P \ 1/2 

bJ j2 \h{Tn,Vr.mn 

\n=Ar+l / 



where Bq = n/2 (^^^^^7=^ j " for q > 2, and = 1 if 9 = 2. 
Taking the expected value on both sides of the inequality, then using 
Minkowski's inequality, we have 



1/9 



E 



J2 Mr„,K)W 



n=N+l 



q/2 



<B^i J2 E[|Mr„,K)Wr] 



912/9 



(3) 

We now, have only to compute E [|/i(r„, y„)(i)|'']. As Vn and r„ are 
independent, as we can compute E[r^'^"] and E[||y„||^] < Mq < 00, 
we have 

T (n- ^ 
E[\h{Tn,Vn){t)\'] < Mq- 



r(n) 



We can use this results in (l3|, leading to 



E[\RN,Pm < B'.Mq [ ^ 

, ii=JV+l 



Tin) 



2/9' 



9/2 



(4) 
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Stirling formulae claims 




r n 



a } 



2/q 



T{n) 



+00 

El a 

n=N+l ^ ' 



which shows us that the series converges. Let us denote 

r in 



~ r(n) 

It can be easily proved that sup„>^^]^ ^n,q — -H^w+i.g, and we have 



(5) 



E 



sup |i?iv.M(i)|' 

Af>JV+l 



< 



with 



a 
2-a 



]S[q/2{2/a-l) ' 
q/2 



□ 



We can improve this result by checking what happens on Ej^ equipped 
with IHI^p. 

Proposition 2.4. For p > and {N + l)a > q > max(p, 2), there 
exists Aq such that 



E 



\Y-YNrK,, 



< 



Vol{Ky'PAq^^^ 



Proof. According to Holder's inequality, we have 

\\Y y^^ll^ p < (vol{Ky-r>/'' (^1^ \Y{t) - dt 



1/p 



where Vol{K) — J^Ik- Thanks to proposition 



2.3 



E 



\Y-Y, 



n\\k. 



< Vol{K) 



q/p-l 



we have 

AqHjM^q 



K 



jV9/2(2/Q-l) 



dt. 



□ 



We have proved that we are close enough to the process Y if we 
truncate the sum. Unfortunately, we are not always able to simulate 
Vn- We will now try to see what happens if we use a sequence of 
random variables {Wn,k)n>i such that limfe^oo Wn.k = in a sense 
we will define later. In the next proposition, we will evaluate the 
distance between Yln=N+i H^n, Vn){t) and J2n=N+i H^n, Wn,k){t) in 

for q > 2. Due to the fact that r„^^" has a q-moment if and only 
if n > q/a, we will not be able to compute in general the distance 
between the sums starting at n 1. 
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Proposition 2.5. For a{P + 1) > a{N +1) > q>2, if {WnM)n>i is 
a sequence of random variables such that there exists a constant Mq^^: 



with /3 < 1/a - 1/2. Then 

1 1 



]Sf2/a-l3-l p2/a-f3-l 



where A'^ depends only on q. 
Proof. Let us denote 



RNAt)= J2 Mr„,K)(i)-Mr„,VF-„,fc)(t). 

n=N+l 



In the same fashion we got Q in the proof of proposition 2.3 
have 



we 



, „=A'+1 



r(n) 



2/9' 



9/2 



Using the same definition of iJ„^g as we did in ([5]), and knowing 
that sup„>jY+i = Hpf^i^q, and using well known series-integral 
comparison we get 



E[\RNAt)n < A'qHN+l^qMq^k 
,9/2 



1 



1 



^2/q-/3-1 p2/Q-/3-l 



9/2 



where A' = ( „ " V 

Q Q \ 2— a/3— a J 



□ 



Proposition 2.6. Forp > anrf (P-t-l)Q; > (A^-|-l)a > q > max(p, 2), 
under the same assumptions as in proposition\2. 5\ we have 



In the same way we had proposition |2.4[ we have 



E 



HN+l,qMq^k 



9 

K,p 



< Vol{K)i/PA' 



1 



1 
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The previous proposition let us evaluate the error we do approxi- 
mating in the sum Ki(t) by Wn,k{t), knowing their L''-distance. In the 
next part, we will use these results by taking enough terms to minimise 
the error done in proposition |2.3| and a good enough approximation of 
Vn such that Mq^k is small, in order to have a small error computing 
J2 h{rn, Wn,k) instead of ^ /i(r„, Vn)- Depending on the values of q 
and a, we may have to deal with the first terms in a particular way. 



3 Application to the local time fractional 
stable motion 

In this section, we will apply the results from section [2] on the process 
defined in the introduction. The precise definition we will work on is: 

where 

• (G„)„>i is a sequence of i.i.d. of standard normal distribution, 

• (X^)„>i is another one, 

• (^n)n>i E^re independent copies of a fractional Brownian motion 
local time, each one defined on some (il', JF', P'). 

• (rn)n>i are the arrival time of a Poisson process of rate 1 on 
[0,oo)7 

all the variables being independent. 

But before that, we will have to prove that this process satisfies the 
assumptions needed. 

In fact, we will be able to generalise a bit this work. Consider 
/„ satisfying the following assumptions: /„ : M x M are the 

independent copies of a non negative continuous random function on 
the probability space {U,'^, T'^, P^) , such that for all if C M compact 
set, denoting Ek the space of continuous functions on K equipped 
with the uniform norm 1 1 • 1 1 , we have for p > 

n\\ux,-)r^] 

uniformly bounded in and fn{-->i) has its support included in 



inf Sf - p, sup B 



H,n 



(7) 



where (i?^'") are independent fractional Brownian motion with Hurst 
parameter H defined on the same space as /„ and p > 0. In the 
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following, we will consider 



(8) 



Remark 3.1. The local time of a fractional Brownian motion satis- 
fies obviously the support condition with p — Q. It satisfies the other 
condition because •) is a non decreasing function so that 



sup ln{x.,t) 
teK 



Inix, to), 



to being the upper bound of K . It simply claims E [||Z„(x, < ^q. 

The following lemma is a direct consequence of the support condi- 
tion so we will skip its proof. 

Lemma 3.2. Let f„{x,-) be a continuous function on K, ||-||^ being 
the uniform norm on K . Let /„ satisfies the assumptions stated above. 
Let Lp be the positive density of Xn distribution with respect to Lebesgue 
measure. For q > 0, if for all a > 0, Jg^ip{x)~'^^"'^^e~'^^ < oo, there 
exists C such that 



E 



^(X„)-i/"/„(X„,.) 



K 



< C sup 



E 



\fn{x,-)\ 



1/p' 



for some p' > 1. 



For some technical reasons that will appear later, we would like 
to work with a slightly different process, having the same distribution 
as Y. Because of that, we will show that the process Y and Y^ are 
identically distributed, where 

Y^{t) = Er-i/"G„^(X„)-'/"/„(^„,i), 

where Xn are i.i.d. random variables, and ip is the density of their dis- 
tribution which satisfies for q > 0, for all a > 0, tp{x)~''/°''^^e~°'^ < 

00. 

In order to achieve that objective, we will show that Y and Y^ have 
the same characteristic function. We can compute this characteristic 
function by using techniques that can be found in |10j . thus we get 

Proposition 3.3. Let K be a compact subset ofR^^, Ex be the space 
of continuous functions on K, equipped with the uniform norm \\-\\, 
and E'x its dual space. For all y' G E'^, we have 



E 



exp / I e 



'<y' ^''^'^'-> -l-i< y',h{r,v) > l|iftM)|l<i) dr\{dv), 
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where A is the distribution ofVn — G„iy9(X„) ^^"/n(-'^n, ■) and 
h X Ek ^ Ek 

Then straightforward calculus lead to the following proposition. 

Proposition 3.4. // for some q > and for all a > we have 
Jg^ip{x)~'^^°'^^e^°'^ < oo, then processes Y^p and Y have the same 
distribution. 

From now on, we will only use 

^(•) = E r-i/"G„e2|^"l/"/„(X„, .), (9) 

n>l 

where X„ has a Laplace distribution of parameters (0,1/2), i.e. its 
density is e"^'^' with respect to the Lebesgue measure on M. 

Proposition 3.5. Let K C denote a compact set. The series 
defining the process Y{t) in ^ converges uniformly on K . 

Proof. The proof of this proposition simply consists in verifying that 
Y satisfies the assumptions of theorem 2.4 in [17 , that is: The measure 
F defined for all borelians in Ek in ([2| is a Levy measure. We have 
to prove 

(i) 



{< y' ,y M)F{dy) < cx) 
(ii) The function 

$ : ^ exp (^j^ (^e'<v'.v> - 1 - z < y', y > l||y||<i) F{dy)^ 

is the characteristic function of a probability measure on Ek ■ 

(ii) is a direct consequence of proposition 3.3 To prove (i), a 
straightforward computation leads to: 

{\\y\\^ M)F{dy) ^ f {v^ X'{dv) (l + 



a 



where the second member of the equality is finite because by indepen- 
dence of the Gn- 

f |z;r A'(d«) =E[|G„nE[(^(X„)-isup/„(X„,i)"]. 

Using lemma 3.2 the integral is finite. Theorem 2.4 in [17] can thus be 
applied. 

□ 
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We will now see how we can apply propositions 2.4 and 2.6 
Remark 3.6. For K C M+ compact, 

Mq :=E[||K|iy <(X3. 

Indeed, by independence of the variables. 



M, = E[|G„r]E 



K 



and the last expectation can be bounded using lemma\37^ 



Using proposition 2.4 we directly get the following corollary 

Corollary 3.7. For K C M+ compact, there is a constant C such that 
for p > and P{a +1) > q > niax(p, 2) 



E 



K,p 



< 



c 



Now, let us study the non-truncated terms. Denoting 
9n,k{^^t)^ / fkiy - x)ln{y,t)dy, 



where {(fk)k is an approximate identity with support in [—\/k,\/k]. 
We will use f{x) = — |a;| + 1 on [—1,1], tp = Q elsewhere. We will 
denote (fikix) = ktp{kx). We can rewrite gn,k as 

gnAx^t)=[ VkiB^'"" - x)ds. 



We will now denote /„ the discretisation of this integral calculated 
using the rectangle method using m„ k points uniformly spread on 

[o,r] 

[m„,fct/T] 



In,k{x,t) 



1 



1=0 



where [x] is the floor function. We have Ki(-) — G„e^l'^"l/"^„(X„, •) 
that will be approximated by Wn^-) — G„e^l"^"l/"/„,fe(X„, •). 
In the following, C will denote a generic constant. 

Proposition 3.8. Let K denote a compact subset o/M+. For q> Q, 

for (3 < \/ a — 1/2, for all 6 < — ^, if we take nin.k — [k^n~y] 
with S' < H, there exists C such that: 

Cln 
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Proof. In this proof C{uj) will denote a generic random variable with 
finite moments of all order. 
According to lemma [3^ 

E[\\Vn-Wn,k\\K] <Csnp(E 
Let us write 

ln{x,t) - In^k{x,t) = {ln{x,t) - gn,k{x,t)) + {g„^k{x,t) ~ In^k{x,t)). 

First, consider of l„{x, t)-gn,k{x, t) = ipkiy-x){ln{x, t)-/„(y, t))dy. 
As the fractional Brownian motion is locally non determinist, we can 
apply theorem 4 in [TS], in order to have for all i5 < jjf ~ 5> there 
exists C{uj) > which has finite moments of all orders such that 

\ln{x,t)-gnAx,t)\<^^. (11) 

Now, let us consider 



\\ln{x, ■) - In.k{x, 



(10) 



gn,k{x,t) - In,k{x,t) 
[m„,fct/T] „jT/m„ 

E 



i=l 



(i-l)T/m„,fc 



MB"''' - x)d.s -ipklB 



We can remark that ipk is fc'^-Lipschitz, and that for S' < H 



ds 



(12) 



sup 

s,teK \t-s\" 

has finite moments of all order. Consequently, using that majoration 
in ( 12 1 claims 



l5n,fe(a;,t) - /„,fc(x,t)| < 



(13) 



6 + 2 _JL-, 

iTT'Ti.k = [k~s^n "5'] with 



Combining ( [ll| and (13) and by taking 
/3 < 1/a — 1/2, there exists C{llj) with finite moments of all order such 
that 



\ln{x,t) - In,k{x,t)\ < 



Using the majoration ( 14 1 in ( 10 1 concludes. 



(14) 

□ 
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Using proposition |2.6| we get 

Corollary 3.9. For p > and {N + l)a > q > max(p, 2), for (3 < 
1/a — 1/2, for all 6 < ^ — |, if we take nin^k — [k^n^'y] with 



S' < H , there exists C such that 
p 



E 



h{T^,Vn)-h{T^,Wn^k) 

n=N+l 



K,p 



< 



C 



1 



kqS ]\Jq(2/a-f}~l)/2- 



For the last part of this section, let us remind some notations 

N 



and Y{t) = limAr^oo YN{t). Let us denote 



N 



(15) 



We will consider ¥{\\Y ~ Zp^kWn p > ^e) for {P+l)a > {N+l)a > 
q > max(p, 2) > p> 1. 



{\\Y ~ ^p,fe||^,p > 3e) < P (||y - yp||^,p > e) 

+ P (ll (yp - y^^) - (Zp,fe - ZN,k)\\K^p > e) 



■(II 



K II A',p 



Corollaries |3. 7| and |3 .9| combined with Markov's inequality allow us 
to evaluate without any difficulties both 



•(II 



y-Yp\\K,r>>^) < 



C 

£gpg(l/a-l/2) 



and 



{\\{Yp - Yn) - {Zp^u - ZnMk^p > < 



C 

eiki^ Ar9(i/a-0-i/2) ' 



Now, we have to study the remaining terms ft.(r„, — /i(r„, Wn^k) 
for n<N. 
Denote 



and 



Cfe(t) = r-i/-G„e2|^"l/"(?„(x, t) - g„,,(x, i)). 
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Since ^ti + Ck = - ^jv,fe(i), we have 



'(II 



Y, 



N 



N 



n=l 



> 



2NJ 



N 

Eip(ll^"' 



> 



2NJ 



(16) 



A direct computation is enough to evaluate 

(2iV)"C 



(lie 



fclL,p > 2iv) - 



For n < N taking rrin^k 



r-p-l/((5'Q), 2+^1 , 

[1 „ ' K i' J we have 



E 



r„ 



q/a 



S'q 



< 



c 



kq(5+2) ' 



SO we can apply Markov's inequality on the last term. ( 13 ) leads to 
the evalution of 



(lie 



Jl >— )< 

k\\K,p 2NJ ~ 



{2NYC 

^qf^qS 



Once this remaks have been made, we only have to tune our pa- 
rameters to control the error before stating the following theorem: 

Tuning procedure 3.10. We want to approximate process Y of pa- 
rameters {H,a) by a family (^p^,fc,)e>o- We see process "K as a series 
and Zp^ k^ as a truncated series. 

We can adjust two parameters: is the size of the truncature and 
kg controls the approximation of the fractional Brownian motion local 
time approximation. We will make two kinds of errors, one coming 
from the truncature itself, and one from the approximation of the local 
time. 

The error from the truncature can be controlled if we take ~ 

The approximation error has two sources: our theorical approxi- 
mation of the local time, and the discretisation used to compute the 
approximation . 

We can deal with the first one by taking fc^ ~ Ce^''/'' where 
5 < 1/{2H) — 1/2 comes from the H5lder continuity of the fractional 
Brownian motion local time. 

Let us denote by m„.fe the number of points used in the discretisa- 
tion. In the series defining Y, the first terms are the most important 
so we will distinguish two cases, the first N terms, where p > 1 and 
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(P + l)a > {N + l)a > q > niax(p, 2), where we will take rrin^k 

-l/(S'a) — 

random, linked with the size of r„, m„,fc ^ [Tn 'h'' ], where 
S' < 1/H comes from the H5lder continuity of the fractional Brownian 
motion. For the remaining terms, a good precision is not as impor- 
tant, so we will need fewer points in our discretisation. We can take 

iTT-n.k = n~s^] with /3 < 1/a — 1/2. 



Theorem 3.11. With tuning procedure 3.1C\ we are able to get a con- 
vergence rate for family of process Zp^^k^ defined in (15). For rj > I 
and p> I there exists C such that 



P(\\Y-Z 



\K,p 



>e)< Ce"^''-!) 



Using P^, and m„ ^ as defined previously, we can have one last 
theorem: 

Theorem 3.12. Zp_,k^ converges almost surely towards Y. 



Proof. For the assumptions used in theorem |3 . 1 1 1 with q > 2 and N > 
q/a — 1 fixed. 

Using Minkowski's inequality: 



E 



1/9 



E\\\iY-YpJ\\^'^' 
-E HYp^_-YN)-{Zp^^k^-ZN,k^ 



Thanks to corollaries |3.7| and |3.9| using the expression of , k^ 



and m„ given in 3.10 we get each term bounded by Ce^, with 77 > 1. 



Thus, we can say that 



E 



\\{Y-Y^)-iZ 



Z 



N 



ks)\\K,p 



1/9 



(17) 



Denote q such as qrj > 1. Inequality (17 1 and Borel-Cantelli lemma 
imply Zp^^k^ — Zpf^k^ converges towards Y — Y^. We only have to 
prove Zpf^k^ converges almost surely towards Y/v. But according to 



(14l, using kg and rUn^k chosen, \\ln{x,t) — In^k{x,t)\\j^ converges al- 
most surely towards so Fat — Z^^k^ too since only a finite number of 
terms is considered. □ 



4 Simulation 

In |3] , the authors explained how to simulate the Brownian motion local 
time stable motion with a random walk with random rewards approach. 
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Figure 1: Allure of Y ioi H = 0.3 and a = 1.2 



They are many advantages of our approach against the random walk 
with random reawards approach. The most obvious is that it is vahd 
for all Hurst parameter and not only H = 1/2. Moreover we already 
highlighted that we have a convergence rate which was not the case 
previously. 

Unfortunately, as one can sec in our tuning procedure, the param- 
eters we use depend on a and H. If a is close to 2, the number of 
terms in the sum is to high to be accepted and if H is close to 1, the 
precision needed in the simulation of the fractional Brownian motion 
is also to high to do. Even if this method is not perfect for all values 
of H and a, it is still a major improvement since we can have sample 
paths for many different values of H with a convergence rate. See for 
example one simulation in figure [T| 

According to theorem 5.1 in [4], the Holder exponent d of our pro- 
cess is such that d < 1 — H. It means that the closer to 1 is H, the 
less regular is our process. See figure [2] for sample paths with different 
values of H and a constant. 

Let us try to make a more rigorous comparison of our two ways to 
simulate our process. We have to take H = 0.5, so that both methods 
are able to do the simulation. In order to check the validity of both 
approaches, we are going to compute E l^e™^^*^] . A straightforward 
computation from the result of proposition [33] gives in the case a — 1, 



E 



,Y{t) 



exp (C \u\ t) 
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Figure 2: Allure of Y for H = 0.2, H = 0.4 and H = 0.6 with a = 0.7 
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Figure 3: Allures of log(E |^e™^*-*^] ) with respect to t. Our method on the 
first line, and the random walk method on the second. 
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We are going to check for both of our processes if log(E [e™^'-*']) is 
a straight Hne using a classical Monte Carlo method. The result is 
shown in figure |3] both methods using 10000 steps, the first one being 
the method developed in this paper, the second the random walk with 
random rewards method. 

According to a linear regression, the E? coefficient is 0.9997 for the 
first one, and B? — 0.9763. Thus we can conclude that even if the 
random walk method seems a bit quicker, our method has not only the 
advantage of being able to deal with H ^ 0.5 but also is closer to what 
we should theoretically get. 

References 

[1] A. Benassi, S. Cohen, and J. Istas. Identification and properties 
of real harmonizable fractional Levy motions. Bernoulli, 8(1):97- 
115, 2002. 

[2] Hermine Bierme, Mark M. Meerschaert, and Hans-Peter Scheffier. 
Operator scaling stable random fields. Stochastic Process. AppL, 
117(3):312-332, 2007. 

[3] S. Cohen, C. Lacaux, and M. Ledoux. A general framework for 
simulation of fractional fields. Stochastic Processes and their Ap- 
plications, 2007. 

[4] S. Cohen and G. Samorodnistky. Random rewards, fractional 
brownian local times and stable self-similar processes. Ann. Appl. 
Probab., 16:1432-1461, 2006. 

[5] ME Crovella and A. Bestavros. Self-similarity in World Wide Web 
traffic: evidence and possiblecauses. Networking, IEEE/ACM 
Transactions on, 5(6):835~846, 1997. 

[6] A. Janicki and P. Kokoszka. Computer investigation of the Rate of 
Convergence of Lepage Type Series to a-Stable Random Variables. 
Statistics, 23(4):365-373, 1992. 

[7] L Kaj and MS Taqqu. Convergence to fractional Brownian motion 
and to the Telecom process: the integral representation approach, 
Brazilian Probability School, 10th anniversary volume, Eds. ME 
Vares, V. Sidoravicius, 2007. 

[8] A.N. Kolmogorov. Wienersche Spiralen und einige andere interes- 
sante Kurven im Hilbertschen Raum. CR (Doklady) Acad. URSS 
(NS), 26(1):15-1, 1940. 

[9] C. Lacaux. Real harmonizable multifractional Levy motions. 
Annales de I'Institut Henri Poincare/Probabilites et statistiques, 
40(3):259-277, 2004. 



17 



[10] C. Lacaiox. Series representation and simulation of multifractional 

Levy motions. Adv. Appl. Prob, 36(1):171-197, 2004. 

[11] M. Lcdoux and M. Talagrand. Probability in Banach Spaces:: 
Isoperimetry and Processes. Springer, 1991. 

[12] J.B. Levy and M.S. Taqqu. Renewal reward processes with heavy- 
tailed intcrrcncwal times and heavy-tailed rewards. Bernoulli, 

6(l):23-44, 2000. 

[13] B.B. Mandelbrot and J.W.V. Ness. Fractional Brownian Motions, 
Fractional Noises and Applications. SIAM Review, 10(4):422-437, 
1968. 

[14] T. Mikosch, S. Resnick, H. Rootzen, and A. Stegeman. Is network 
traffic approximated by stable Levy motion or fractional Brownian 
motion. Ann. Appl. Prohab, 12(l):23-68, 2002. 

[15] L. Pitt. Local times for gaussian vector fields. Indiana Univ. 
Math. J., 27:309-330, 1978. 

[16] J. Rosinski. On Path Properties of Certain Infinitely Divisible 
Processes. 1987. 

[17] J. Rosinski. On Series Representations of Infinitely Divisible Ran- 
dom Vectors. The Annals of Probability, 18(1):405 430, 1990. 

[18] J. Rosinski. Series representations of Levy processes from the 
perspective of point processes. In Levy processes, pages 401-415. 
Birkhauser Boston, Boston, MA, 2001. 

[19] G. Samorodnitsky and M. Taqqu. Stable Non-Gaussian Random 
Processes. Chapman & Hall, 1994. 



18 



